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Abstract. The theory of affine processes on the space of positive scmidcfinitc 
dxd matrices has been established in a joint work with Cuchicro, Filipovic and 
Teichmann (2011). We confirm the conjecture stated therein that in dimension 
d > 1 this process class does not exhibit jumps of infinite total variation. This 
constitutes a geometric phenomenon which is in contrast to the situation on the 
positive real line (Kawazu and Watanabe, 1974). As an application we prove 
that the exponentially affine property of the Laplace transform carries over to 
the Fourier-Laplace transform if the diffusion coefficient is zero or invertible. 



1. Introduction 

Affine processes are a special class of stochastically continuous Markov processes 
with the following feature: Some suitable integral transform (such as the charac- 
teristic function [5], Laplace transform ([EH3]), Fourier-Laplace transform, or even 
moment generating function [7]) of their transition function is exponentially affine 
in the state variable. It has become custom to describe affine processes in terms 
of a parametrization of their infinitesimal generator quite similarly to the Levy 
class [14], where the so-called Levy-Khintchinc triplet (a, c,m(<i£)jj relative to a 
truncation function x(£) allows a parametric description of the generator 

Af(x)=af"(x) + bf'(x)+ J (f(x + 0-f(x)-f(x)x(0)m(<%). 

K\{0} 

The affine property translates into affine drift, diffusive and jump behavior, and 
the coefficients of the involved affine functions determine the so-called "admissible 
parameter set"([S]). For instance, for the state space R+ := [0, oo), Kawazu and 
Watanabe (8j show that the infinitesimal generator of a conservative affine processes 
X takes the form0 

Af(x) = axf"(x)+(b+l3x)f'(x)+ J (/(^+C)-/(x)-/'(x) x (0)(™(^)+^(^)), 

R+\{0} 
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with "parameters" (a > 0, b > 0,/3 e E, m(d£), where the last two objects 

are sigma-finite measures on R+ \ {0} such that 

Ai)m(de)<oo, / (iief ai)m^)<oo. 

^R+\{0} 
K+\{0} 

However, a full parametric characterization depends crucially on the geometry of the 
state space, and the probabilistic properties of afhnc processes may vary accordingly. 
Motivated by multivariate extensions in the afhnc term structure literature as well 
as in stochastic volatility, Duffic, Filipovic and Schachermayer [5] establish a unified 
theory on the so-called canonical state spaces M" 1 x R ra (for further insights, and 
certain simplifications, see OUT]). The recent theory of Cuchiero, M., Filipovic and 
Teichmann [JJ for afhnc processes on positive semidefinite matrices is a response 
to suggestions in the finance literature concerning affine multi-asset models based 
on matrix factors. Those, in turn, have mostly used the class of Wishart processes 
as put forward in [3], or the OU-type processes driven by matrix variate Levy 
subordinators [2]. For a review on financial modeling issues with matrix factors, 
see the extensive introduction of [3], as well as the references given therein. 

Aim of this paper is to show that affine processes on Sf, d > 2, do not exhibit 
jumps of infinite total variation (Theorem I3.2[) . This important result confirms a 
conjecture formulated in [JJ Section 2.1.4]. In the conservative case, it allows to 
simplify the semimartingale decomposition of [H Theorem 2.6]; this is subject of 
Theorem 13.41 A crucial application of Theorem 13.21 concerns the affine character of 
the Fourier-Laplace transform of affine processes (Theorem 14. ip . In particular, we 
show that in the presence of non-degenerate diffusion components, affine processes 
are affine in the sense of Duffie, Filipovic and Schachermayer (Theorem 14.21 and 
Corollary |4.3[) . This means, that their characteristic function is exponentially affine 
in the state- variable. A detailed introduction to this topic with technical remarks 
is given in Section |U 

The main result of this paper, Theorem 14.11 reveals a geometric phenomenon; 
Indeed, in the much simpler case d = 1, where the state space simplifies to the 
positive real line R + , stochastic processes with jumps of infinite total variation 
actually exist. For instance, let <i£ denote the Lebesgue measure on [0, 1] and 
define a linear jump characteristic as 

p(dO :=r 2 l(o,i](0# 

Clearly 1 = f Q £, 2 n(d£) < oo, hence due to [5], an affine pure jump process X with 
infinitesimal generator 



Af(x):=x J (f(x + S)-f{x)-f\x)Z)n(dS) 
R+\{o} 

exists. Let us start X at x > and denote by X t its cadlag representation^. Then 
X t is a special semimartingale with characteristics (A = 0, B = 0, v{dt, d£)), where 
the compensator of Xt equals 

v{dt,d£) = X t n{d^)dt. 



^Such exists due to the Feller property of X, see \E\. 
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The canonical decomposition of Xt is given in terms of the Poisson random measure 
associated with its jumps, /j, x (dt, d£): 

i 



X t = x + £* (p A - v) = lim / / £(/i A {da, d£) - v(da, d£)) 

?>£ 

and clearly X t > a.s., because the jumps of X are positive throughout. Hence, 
almost surely it holds thaO 



]T \AX S = ^2 Ax s = J £ KdO x (J X s ds 

s<t s<t R+\{0} 



For d > 2 the complex geometry of the boundary dS^ of Sf ~ it is not anymore 
the origin only, nor it is a smooth manifold - leads to non-trivial restrictions con- 
cerning the linear jump behavior. One of these is (|3.3j) below, which expresses that 
transversal to dS^ only finite variation jumps are allowed. In addition, there is a 
non-trivial tradeoff between linear drift and linear jumps, see eq. (|3.6|) . One of the 
nice consequences of Theorem [321 is that these two conditions may be disentangled 
from each other, into a simple condition that the drift must be inward pointing 
at the boundary (eq. (|3.4[) ) and the compensator of the affine processes satisfies a 
stronger integrability condition (see (|3.3[1 ). Furthermore, the admissible parameter 
set is now formulated independent to truncation functions, which is impossible in 
the setting of canonical state spaces [5], and in particular for d = 1. 

It should perhaps be noted that the original characterization of affine processes 
[H Theorem 2.4] and all consequences thereof arc stated in a way, which nest 
the one-dimensional one (cf. [8] and [5]) - this is possible in view of the implicit 
nature of condition (|3.6[) . As such the preceding theory is perfectly valid in its 
original formulation; the contribution of the present work, however, is a technical 
simplification of the theory of affine processes on the cone of positive semidefinite 
matrices of arbitrary dimension d > 2, and the additional theoretic results 
concerning the Fourier-Laplace transform of this process class. 

2. Notation and definition of the affine property 

We try to keep notation and presentation of this paper as simple as possible. 
As reference, both for applied and theoretic issues, see the quite extensive work 
[4]; this also concerns technically involved facts, which are here only recollected in 
prose. 

1a equals the indicator function corresponding to some set A. Sd denotes the 
linear space of real d x d symmetric matrices, and (x, y) := tv(xy) is the standard 
scalar product thereon, given by the trace of the matrix product. Accordingly, || 
is the induced norm on Sd, and the pierced unit ball equals 

Bl := {z eS+\0< \\z\\ < I}. 

The natural order introduced by the closed convex cone Sj" is denoted by ^. The 
cone of positive definite matrices is denoted by S^ + (and clearly is the interior 
of S^). V/(x) denotes the Frechet derivative of a function / at x E Sd- This 
coordinate free notation allows much a much shorter and more elegant presentation; 



4 Of course in the finite case, the two summands would differ in general. 
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for an account of the involved details and the coordinate wise way to write what 
follows, the reader is referred to the nice exhibition p] as well as [4]. Only in the 
proof of the main theorem 13.21 coordinates are used. 

We consider a time-homogeneous Markov process X with state space Sf and 
semigroup (Pt)t>o acting on bounded Borcl measurable functions / 



Here pt {x, d£) denotes the (possibly sub-)Markovian transition function of X. 

Definition 2.1. The Markov process X is called afhne if 

(i) it is stochastically continuous, that is, lim s _>t p s (x, •) = Pt{x, •) weakly on 
Sj" for every t and x G S^, and 

(ii) its Laplace transform has exponential- affine dependence on the initial state: 



for all t and u,x G Sj , for some functions <f> : R + x S'J — > R+ and 



The so-called admissible parameter set is introduced in the following. Note that 
unlike [H Definition 2.3] truncation functions may be omitted, and the complicated 
admissibility condition (|U (2.11)], see also (|3.6|) in the proof below) is dropped: 

Definition 3.1. Letd > 2. An admissible parameter set (a,b,B,c, r y,m(d£),fj,(d£)) 
consists of 

• a linear diffusion coefficient a G Sj", 

• a constant drift term b G Sf which satisfies 





(2.1) 



3. Main Result and Proof 



bh{d- l)a, 



(3.1) 



• a constant killing rate term c G R+, 

• a linear killing rate coefficient 7 G , 

• a constant jump term: a Borel measure m on Sj" \ {0} satisfying 




(3.2) 



si\{o} 
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• a linear jump coefficient [i which is an -valued, sigma finite measure 



S+ \ {0} andE 

A 1) < oo. (3.3) 



on 



sj\{0} 

• and, finally, a linear drift B, which is a linear map from Sd to Sd and 
"inward pointing " at the boundary. That is, 

(B(x), u) > for all x,u £ Sj with (x, u) = 0. (3.4) 

The main statement of this paper follows: 

Theorem 3.2. Let X be an affine processes on S^ (d > 2). Then its Markovian 
semigroup (Pt)t has an infinitesimal generator A acting on the space of rapidly 
decreasing function^ supported on St 

Af(x) = 2<VaV/(z), x) + (Vf(x),b + B(x)) - (c+ ( 7 , x))f(x) 
(f(x + 0- f(x))(m(dO + x)). 

si\{o} 

where (a, b, B, c, 7, to, /z) are an admissible parameter set in the sense of Definition 
\3J\ 



Proof. Let x(£) : S^ —> S'j" be a truncation function, that is x(0 = ? near zero, and 
X is continuous, and (what may be assumed without loss of generality) bounded by 
1. By [H Theorem 2.4 and Definition 2.3] the semigroup (Pt)t has an infinitesimal 
generator A acting as 

Af(x) = 2<VaV/(z), x) + (Vf(x),b + B(x)) 

-(c+( 7 ,x))f(x)+ J (f(x + Z)-f(x))m(dO 

S d +\{0} 

+ I (/(i + 0-/W-MO,v/(i)))W<i£) 1 i). 

S d +\{0} 

where c > 0, 7 £ Sj~ and the parameters a, b £ S^ satisfy (|3.1[) and is a sigma 

finite S'j" valued measure on S'j"\{0} which integrates ||£|| 2 A1. Furthermore, m(d£) 
is a Borel measure on 5'j"\{0} which satisfies (|3.2[) . We also know from [4j Theorem 
2.4 and Definition 2.3] that jumps of infinite total variation may only occur parallel 



^We deviate here a little from the corresponding definition in [4], where /i is a finite measure 
on Sj~ \ {0} (later divided by ||§|| 2 A 1). Here "S^-valued" has to be understood as follows: For 
any Borel set E in £j~ such that its S^-topological closure E C \ {0} we have n(E) £ . 
This allows infinite activity jumps with state-dependent intensity: Indeed, there exist fi for which 

\{0}) = 00, which nevertheless satisfy eq. 1 13.311 . The latter simply means that (||£|| A 
is a finite measure. 

^The integral is of course matrix valued, and < 00 expresses that it is finite. 
^For further details, see [4] Appendix B]. 
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to the boundary. In terms of admissibility conditions of the parameters, this is 
expressed in eq. ()3.2[) as well as the following condition: 

J {£,u}{fi(d{;),x) < co, x,u G with (x, u) = 0. (3-5) 

S d + \{0} 

Furthermore, the drift must be inward pointing at the boundary. That is expressed 
in the terms of the positivity of the constant drift b 6 Sf (as above) , as well as the 
following constraint on the linear part B (a linear map from Sd to Sd)' 

(8{x), u) - J (x(0, «> {KdO, x}>0, x,ue S+ with (x, u) = 0. (3.6) 

Note that (x,u) — is equivalent to xu = ux = 0, that is x, u G dSt, see also [4j 
Lemma 4.1 (i)-(iii)]. 

Suppose for a moment that the validity of (|3.3|) had already been shown. Then 
\x integrates ||x(£)ll ^ ll£ll A 1 an d therefore a new drift may be introduced as 

B(x) := B{x) — J X (0(KdO,x) (3-7) 
st\{o] 

which, in view of ()3.6[) satisfies admissibility condition p.4p . Hence the proof of 
the theorem were settled. So the essential point of the statement is (|3.3[) . We use 
standard Euclidean coordinates on Sd for the remainder of the proof; all indices 
range between 1 and d, if not otherwise stated. Let {c l - J , i < j} denote the canonical 
basis of Sd, that is, the kith coefficient of c y equals 

c l 3 i = hkSji + SjkSu(l - Sij), 

where Sij denotes the Kronecker delta. If i = j then = SikSn (the diagonal 
matrix with a 1 in the ith diagonal entry and zeros everywhere else) . Otherwise c lJ 
is zero except in the zjth and jith entry, where it is equal 1. We may evaluate /i 
coordinate wise as fx{A) = (fiij(A))ij, A G Sd such that ^ A, the latter denoting 
the topological closure of the set A. Let c* := 1 — c M , where / is the unit matrix. 
Then clearly (c u , c|) = 0, and by eq. ((33]) it holds that 

| C/'u (dO < oo, 1 '../ • '/• ' •' j. (3.8) 

Bf 

We show now that a similar intcgrability condition must also hold for 1 < i = j < d. 
To circumvent intcgrability issues at the origin, the measure \i is pierced as follows 
near 0: For e > we introduce the new (and by construction finite measures) 
[f{di) := Ai(^) 1 e<||?||<i( rf ?)- In particular, 

Vij{d£) = ^l e <||cj|<i(^) 
arc all signed finite measures (1 < i,j < d) such that for all e > we have 

- oo < / Cfe/MljK) < oo, 1 < i, j, k,l<d. (3.9) 
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By (|3.8[) . there exists a positive constant M such that for all e > 



0< / CulijM)<M, i^j. (3.10) 



"33 

S+ 



We introduce now the following boundary points of St. 

e% := c u ± c ij + 6> j , i^j 

By construction (eY,ej?) = 0. Setting x = eY and u = e°_ and applying (|3.5p . we 
must have 

< J (& - 2& + + 2/iij(df) + Mii(df)) < oo, * / i- 

Similarly, we obtain by using x = e!f and u = e+ that 



< / (& + 2& + & ) (M« (df) - 2/iij- (df ) + Mij (df )) < oo, i ^ j. 

Accordingly there exists a constant positive Mi we have, for all e > 0, 

< J (& - 2& + & )(/£(#) + 2/4 (d£) + < Mi, i ± j. (3.11) 

and 

< j (& + 2& + &j)04(#) - 2^-(df) + < M x , i ^ j. (3.12) 
Summing up (|3.1ip and (|3 . 12[) we therefore obtain 

o< / {^pUdO-^i^U^ + ^j^O) 



+ / (f«Mji(de)-26i/4j(df)+6i/*«(«)) <2Mi 



for all i ^ j. The two integrals arc non- negative, because ^ e is an S'j" valued 
measure. We therefore conclude that both of them arc finite: 



< J (Siil4i(d0 - 2^-4-(d0 + tnl'jA'lO) < 2Mi, i ^ j (3.13) 

o< / (e«/i|i(de)-2f«M«(de)+M(de)) <2Mi, (3.w) 



By subtracting (|3.10j) from (|3.14j) twice, once for i, j and then j, i, we have 
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for all e > 0. Plugging this information back into ()3.13j) and using the fact that 
£a > 0, and \x e is positive semidefinite, we obtain 

< / (f«M«(dO + < 2(Mi + M). 



The choices of z was arbitrary. Taking into account ()3.8[) and the preceding uniform 
estimate in e, we finally conclude 



< / i„ /' , ,•:'/: ) < oo, l<i,j<d. (3.15) 

Define the positive measure tr(/i)(d£) on Borel sets A with ^ A by tr(fi(A)). 
Eq. (|3.15|) implies immediately 

tr(£)tr0u)(O < oo. (3.16) 



We finally show the admissibility condition f|3 . 3[) : Let £ be a positive semidefi- 
nite matrix with diagonalization £ = UDU T , where f7 is orthogonal and I? = 
diag(Ax, . . . , Ad). By using this diagonalization and the cyclic property of the trace, 
we obtain 



d 



\C\\ 2 = tr(UD 2 U) = tr(D 2 ) = £ A 2 < A,) 2 (3.17) 

i=l i=l 

= (tr£>) 2 = (trUDU) 2 = tr(£) 2 , 



where < follows from the non- negativity of the eigenvalues Ai. Using this technical 
detail, we infer from (|3 . 1 6[) the following estimate: 

f (Hell A l)tr( M )(dO < oo. (3.18) 
Js+\{0} 

By Lemma 13.31 below we may conclude the validity of condition (|3.3|) . Hence the 
definition of B by cq. (|3.7p is legitimate. □ 



The following technical statement has just been used and again will be used in 
the proof of Theorem 14.11 

Lemma 3.3. For any non-negative B or el-measurable function g we have 

II / gitMdOW < [ 0(Otr(/*)(de)<oo. 
Js+\{o} Js+\{o} 

Proof. Since \x is a positive semidefinite measure, we have by eq. (|3 . 1 T[) for any 
Borel set B C S% \ {0} the estimate \\n{B)\\ < ti(n(B)). Hence approximating the 
function g by non-negative simple functions, the assertion follows. In particular, 
the integral (|3.3[) must be finite whenever (|3 . 1 8|) is finite. □ 



AFFINE PROPERTY AND JUMP BEHAVIOR 



9 



3.1. The Semimartingale Decomposition. Suppose X is a conservative affine 
process on Sj", d > 2, with admissible parameter set (a, b, B, 0, 0, m(d£), fi(d£)). In 
view of the Feller property [H Theorem 2.4] of X, for each initial state x € 
there exists a modification of X t := (X x ) t >o on the canonical path space, which 
is a cadlag semimartingale. Since we know that the jumps of X are of total finite 
variation, we have as an immediate consequence of [H Theorem 2.6], 

Theorem 3.4. Let S be a d x d matrix such that S T S = a. Then there exists, 
possibly on an enlargement of the probability space, a d x d-matrix W of standard 
Brownian motions such that X admits the following representation 
t t t 

X t =x + bt + J B(X s )ds + J (^/Y s dW s T, + T^ dWs^Cj + J J £ii x (ds,d£), 
oo S d 

where pL X (d£,ds) is the random measure associated with the jumps of X , having 
compensator 

v{dt,d£) = (m(df) + (X t ,fi(di)))dt. 

Note that if the drift is of the particular form B(x) = fix + x(3 T , where j3 is a 
real d x d matrix, if b — 8a (5 > 0) and in the absence of jumps, X is a Wishart 
process ( [fl ES ) . 

4. The Fourier-Laplace transform of affine processes 

Affine processes on positive semidcfinite matrices are defined in terms of the 
Laplace transform of their transition probabilities, cq. (|2.1[) . In general, the Laplace 
transform is a natural choice for integral transform of generalized functions on 
proper cones such as Sf . However, Duffie, Filipovic and Schachermayer [5] have 
defined affine processes on Rip x K™ in terms of the exponentially affine form of their 
characteristic function. Only in the one dimensional case R + , the two state-spaces 
coincide and therefore also the two definitions of the affine property, either by the 
Laplace transform (Kawazu and Watanabe [5]) or by the characteristic function. 

Therefore, the question whether the characteristic functions of a positive semi- 
definite affine process is indeed exponentially affine in the state, is of considerable 
interest. We will denote this property as being "affine in the sense of Duffie, Fil- 
ipovic and Schachermayer" . 

Unless the diffusion coefficient a vanishes, X need not be infinitely divisible, 
or equivalently, infinitely decomposable (for the definition and characterization of 
these properties in the affine Markov setting, see [U Definition 2.7, Example 2.8 
and Theorem 2.9]). This complicates the problem of extending the affine formula 
cq. (|2.1j) to the complex domain, because it is not anymore guaranteed that the 
Fourier-Laplace transform of X exhibits no zeros, as is in the infinite divisible case 
( |14l Theorem 25.17]) (which is a necessary condition to write it in an exponentially 
affine way). From the ODE perspective, there is a related technical problem, namely 
to show that the real part of ip(t, u) as solution of the system of generalized Riccati 
equations (eqs. (|4.1[) - (|4.2[) below) with imaginary initial data does not explode in 
finite time. Indeed, we have the estimate 

| e -<£(t,u)-(V>(i,«),x)| < e -R e (^(t 1 «))-(R e (^(t,-u)),x) 



^For sufficient and necessary condition of conservativeness, see [4] Remark 2.5] and 1121 Section 

3]. 
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and if the real part of ip explodes in finite time, then the characteristic function 
must have a zerod. In this section we extend the affine transform formula to the full 
Fourier-Laplace transform, under the premise that the diffusion component must 
be non-degenerate or equals zero. For technical difficulties in the degenerate case, 
see Remark 14.41 

We denote by <S(Sj~) the complex tube S^+iSd, and similarly S(Sj + ) = + 
iS d and <S(M+) = R + + iR. 

Theorem 4.1. Let X be an affine process on S'j" (d > 2), with a diffusion coefficient 
a which is either invertible or zero. Then the affine property ((2.1)) holds for all t > 0, 
x G Sf, and for all u G S(S^ + ), with exponents <f> : R+ x S(S^ + ) — > S(R+) and 
ijj : R_|_ xS(S^ + ) — > S(S^ + ) which are the unique global solutions of the generalized 
Riccati differential equations 

d t <f>(t,u) = (b,ijj(t,u))+c- J ( e -W*.«).0 - 1) m (d£), (4.1) 

S d + \{0} 

dtip(t,u) = -2ip(t,u)ail>(t,u) + B T (tp(t,u)) + 7 (4.2) 
- | ( e -<*(*.«).€> - l) 
sjAfo) 

given initial data <p(Q, u) = 0, ^(O, u) = u. 

For the following two results we assume as in Theorem 14. 11 that d > 2, and the 
diffusion coefficient a of X is either invertible or zero. 

Theorem 4.2. The affine property (|2.ip ofeo ZioWs /or u G 5(S'j"), exponents 
cj> : R + x — ► 5(K+) and V- 1 : x — » S(S^), being (not necessarily 
unique ) solutions of the generalized Riccati differential equations (|4.ip - (|4.2p . 

Applying the above to u G iSd, we finally obtain: 

Corollary 4.3. X is affine in the sense of Duffie, Filipovic and Schachermayer 

4.1. Proof of Theorem Hm 

Proof. Let u = v + iw G S(S^ + ). We denote by tp(t,v) the unique global solution 
of equation (|4.2p on which exists due to [U Proposition 5.3]. tp(t, u) is defined 
as maximal solution of (|4.2p on the open domain 5(S'j _+ ). Note that the right side 
R(ip) of eq. (|4.2p is an analytic function thereon, hence it is in particular locally 
Lipschitz. Accordingly, the maximal life time of ip(t,u) equals 

t + (u) := lim in£{t > | Reip(t,u) G dSf or ||^(*,«)|| > n} 

n— »oo 

and we have < t + (u) < 00. 

First, we show that ip(t,u) does not touch the boundary of S(S^ + ) in finite 
time. To this end, we introduce the function u) := Re(^>(t, u + iv)), which is 
well defined for t G [0,t + (u)) and has values in . Denote by u i-> R{u) the 



We note that the tf> coefficient does not matter here: Re(</>(t, u)) > can be inferred from the 
specific form of the generalized Riccati differential equations. 
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function on the right side of (|4.2j) . By straightforward inspection, one observes 
that for all t < t+(u) 

d tX (t, u) - Be(R(x(t, «))) t = 9trl>(t, v) - R^(t, «)). 

Since R is an analytic and quasi-monotone increasing function on 5j~ + with respect 
to the cone S'j" (see [U Definition 4.7 and Lemma 5.1], we may invoke Volkmann's 
comparison result in the fashion of [4j Theorem 4.8] and derive 

X(t,u) >z ip(t,v), for all t < t + (u). 

But tp(t,v) e Sj + , for all t > 0, pfl Proposition 5.3]. Hence we have shown that 
ip(t, u) does not touch the boundary of S(S^ + ), which is dS^ x iSa, in finite time, 
and therefore we have 

t + (u) := lim inf{i > | U(t,u)\\ > n}. (4.3) 

n—>oo 

Hence it remains to show that ip(t, u) does not explode in finite time. Since affinc 
transformations of the state space do not effect the blow-up property, we may 
without loss of generality assume that the diffusion coefficient equals zero or equals 
the identity matrix. To obtain the necessary transformation, one can adapt [4j 
Propositions 4.13 and 4.14]. We introduce the shorthand notation K := K\ + K2, 
where ^ 

0<||£||<i 

and 

K 2 (u):= J (l-e-^)n(dO 
IISII>1 

Using this decomposition, we can write 

R(u) = -2uau + B T u + 7 + K(u). 

Using the Cauchy-Schwarz inequality, Lemma 13.31 and condition (|3.3[) , we infer the 
existence of a constant C\ > such that for all u £ S(S^) 

\\Ki(u)\\ < \\u\\ [ (U\\ A l)tr( M )(dO - dH. (4.4) 
^+\{o} 

The same condition allows to conclude the existence of a positive constant C2 such 
that 

\\K 2 (u)\\ < f 2tr(j0(d£) = C 2 < 00 (4.5) 
J||?ll>i 

where we once again have used Lemma 13.31 
By Lemma |B. II in Appendix [B] we have that 

Re(^(t, u), ip(t, u)aip(t, u)} > (4.6) 

for all t < t + (u). Using estimates (14.4[) - (14.6[) and the Cauchy-Schwarz inequality, 
the existence of a positive constant C follows, such that for all u S S(S^ + ) and 
t < t+(u), 

d t U{t, u)\\ 2 =2 Re@(t, u),R(i>(t, «)> 

< 2Re(^(t,u),B T (iP{t,u)) +j + K(i/)(t,u))) 
<2C(l + \\4it,u)\\ 2 ). 
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Hence, by Gronwall's Lemma (or, equivalently, by standard comparison for scalar- 
valued ODEs) we obtain for all t < t+(u), 

\m,u)\\<e ct ^TTW (4-7) 

which in view of (|4.3|) proves that t+(u) = oo. So we have shown that t i— > ip(t,u) 
is the global solution of (|4~2j) for all u G <S(5j"). Moreover, Re(ip(t,u)) G for 
all t > and the right side of (|4.1[) is well defined for all u G S(St). Therefore 
plugging ip(t,u) into (|4.1[) and integrating with respect to time yields <fi(t,u). 
Now for each t > 0,x G $1", the Fourier-Laplace transform 

g(u) = E[e- {u ' Xt) \X = x] 



and the function 

f(u) := g-^f*.")-^^^'")) 

are complex analytic functions on <S(S't + ), and (in view of (|2.1|l ) they coincide 
on set of uniqueness, namely S'j"" 1 "- Hence f = g on. S(Sj[ + ), which proves the 
assertion. □ 



4.2. Proof of Theorem 14.21 We can write u = v + iw, where v G Sj" and denote 
for each n > 1 the matrix u n := (u + ^1) + iw, where 1 is the unit d x d matrix. We 
further denote by i(> n (t) the solution of (|4.2[) . subject to il> n (0) = in, which exists 
globally due to Theorem 14.11 because now u n G S(S^ + ). 

Let 7r(x) be the projection of x G onto Sj~, which exists uniquely, because 
is a closed convex set. For u = v + iw G S(Sd), we slightly abuse notation and 
write 

ir(u) := 7r(t>) + iu> G S(S^) 

Using the continuity of the right sides of (|4.ip ~ (|4.2p we may as well consider <f> n {t) 
and ipn(t) as solutions to the generalized Riccati differential equations 

d t <j> n (t,u) = (b,i> n {t,u))+c- J ( e -«iMt.»O).0 - l) m (de) (4.8) 

■S d +\{0} 

9tV„(t, u) = -2^„(t, u)a^ n (i, tt) + B T {^ n {t, «)) + 7 (4.9) 

- y ( e -H^^»'«) - 1) M (d0 
s+\{o} 

subject to Vn(O) = ««, </>n(0) = on the whole domain S(Sd)- 

Now by estimating (|4.7|) in the proof of Theorem 14. 1[ there exists a uniform 
constant C, such that for each n, 

HVnWII <e c Vl + KII 2 . (4.10) 
But this means that for any T > 0, the family of curves 

{Mt) \t G [0,T]} 

lie in a single compact set if. Since T is arbitrary, an application of Lemma IA.2I 
therefore yields that there exist functions t h-> <p(t,u),t *-> 4>{t,u) on [0, oo) which 
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are the pointwise limits of a subsequence (4>n k (t), ipn k (t)) (k — > oo) and they satisfy 
eqs. (|4.1j) - ()4.2|) . Furthermore, we have by dominated convergence, 

fc— >oo A;— >oo 

= E[e-<"' Xt) | X =x}. 

This ends the proof. 

Remark 4.4. • It can easily be seen either by numerical experiments or 

explicit calculations that (an appropriate adaption of) Lemma IB. II docs 
not hold, if a not equals a scalar multiple of the unit matrix. To be more 
precise, in general, the real part of 

tv(xxax) = tr(xxxa) 

can be strictly negative. For instance, using 





) 




r 1 i \ 


i 


) ' 




, i 4 ) 



we obtain Re(x) £ S£ , but Hetr(xxax) = — 1 < 0. As a consequence, 
we cannot derive estimate (|4.6[) . which in turn is a technical necessity to 
obtain the a-priori estimate (|4.7j) rcsp. (|4.10j) . However, we conjecture 
that the problem concerning degenerate, nonzero diffusion coefficient a 
admits the same answer as that of Theorem 14.21 

• It should be noted, that the main technical complication concerning jumps 
in prior research had been the presence of a truncation function in the 
right side of eq. (|4.2[) . Only the finding of Theorem l3.2l lct us establish the 
general a-priori estimate (|4.7p resp. (|4.10p . 

4.3. Examples with degenerate, nonzero diffusion. In the presence of a non- 
zero diffusion component a, Theorem 14.11 requires that a is invertible. It should, 
however, be reported that if X is "Wishart with state independent jump behavior" , 
then not only is it evident that X is affine in the sense of [S], but also the affinc 
character of the Laplace transform can be extended to the domain S(S^). And 
in this case, we can solve the Riccati equations explicitly, with no non-degeneracy 
assumption on a. 

Definition 4.5. A matrix-variate basic affine jump-diffusion X on Sj" (MBAJD 
in short) is an affine process with parameters 7 = 0, c = 0, fj, = 0, a constant drift 

, „ d- 1 

b = 2pa, p > 

and with a linear drift B of the particular form 

B(x) = f3x + x(3 T , 

where /? is a real d x d matrix. 

Remark 4.6. • If d = 1, and m(dt;) is a multiple of the density of an 

exponential distribution, then X is a BAJD as introduced by Duffle and 
Garleanu in [6] . 

• If d > 2 and m = 0, then X is a Wishart process, see pQ, [3] and [13j . 
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It is quite straightforward to check that any MBAJD is a conservative Markov 
process ([H Remark 2.5]) and that eqs. (|4.1|) - (|4.2|) take the particular form 

d t( />(t,u) = 2p(a,iP(t,u))- J (e-W*^- c > - l) m(de), 

s+\{o} 

d t ip(t, u) = -2ip(t, u)onp(t, u) + ip(t, u)j3 + /3 T ip(t, u). 
In the following we denote by wf the flow of the vector field fix + xf3 T , that is, 

u" : R X 5+ ^ S+, 4{x) := e^xe^K 
Its twofold integral <rf : — > for t > is denoted by 

ct' 3 : l + xS+4 5+, of (a) = 2 / wf (x)da. 

Jo 

By matrix analysis, we obtain the following semi-explicit solutions for initial data 
u e 5(5+), 

0(t,«) =plogdet (/ + uof (a)) - / (c - ^^> c > - l) m(df), 

S d +\{0} 

Appendix A. Convergence of ordinary differential equations 

The following results are consequences of standard ODE theory. The first one is 
clearly elaborated in [TU1 Lemma 8] , and the second one is a variant of [TU1 Lemma 
9] and can be proved similarly as in [TU] (the difference being that we drop the 
Lipschitz continuity of /, hence one cannot show that every involved subsequence 
in IA.2I converges, let alone to the same limit). 

We recall them here for the convenience of the reader, and without any proof. 
We consider the system of ordinary differential equations on K m , 

d t m = f(t,m), (a.i) 

subject to an initial condition ^(0) = u 6 M m . Recall that equation (|A.I[) possesses 
a maximal solution on a half open interval [0,t+(u)) if the function f: I x U — > E 
is continuous. Note however that such a solution may be not unique if / is not 
locally Lipschitz continuous. 

Lemma A.I. Let U C E be open. Let /, /i,/2, • • • be continuous maps from I x U 
to E. Suppose f is locally Lipschitz and f n converge to f uniformly on all compact 
subsets of L x U . Let ip n <E C 1 ([0,6'„), U) be maximal solutions of 

d t ipn(t) =/„M„(*)) (A.2) 

such that ip n (0) converge to some u G U as n — > oo. Then we have 

t+(u) < lim0„. (A.3) 

Let < a < t + (u) and uq be such that 9 n > a for n > hq. Then the sequence 
il>n +k(t), k = 1, 2, . . converges to ijj(t) uniformly on [0, a] as k — > oo. 
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Lemma A. 2. Let U C R m be open. Let /, /i, /2, • • • be continuous maps from 
I xU to R m . Suppose /„ converge to f uniformly on compact subsets of L xU. Let 
< a < T and ip n £ C 1 ([0,a], U) be solutions of (|A2[) such that ip n (0) converge to 
some u £ U as n — > oo. If for some compact set K C U, tp n (t) £ K for all t £ [0, a], 
then there exists a (not necessarily unique) solution ip(t) of equation \A . 1\) on [0,a], 
and a subsequence tp nk (t) — > ip(t) and dti() nk (t) — > dtipif) uniformly on [0,a]. 

Appendix B. A simple matrix inequality 

Lemma B.l. For any complex valued m x n matrix a and for any b £ S(S^) we 
have that 

Retr(6a T a) > 

Proof. Write a = a\ + za 2 , and b = b\ + ib-2- Then we have 

Retr(6a T a) = Retr((6i + ib-i)(a[ — iaj )(ai + ia^)) 

= Retr((6i + ib2)(aj a± + iaj 02 — iaj a± + 0^0,2)) 

= ti(bi(aj ai)) + tr(6i(a 2 02)) + tr(6 2 ajai) — tr (b 2 aj 02) 

> + + tr(6 2 02 °i) _ tr(6 2 r a7a 2 ) = 0. 

Here we have used that ajai, 61 £ 6 2 = 6^ and the commutativity of the 
matrix trace. □ 
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